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ABSTRACT 

Several approaches have been used to search for dark matter in our galactic 
disk, but with mixed results: maybe yes and maybe no. The prevailing ap- 
proach, integrating the Poisson-Boltzmann equation for tracer st ars, has led to 



more definitive results: yes and no. The touchstone yes analysis of iBahcall et al. 



(119921 ) has subsequently been confirmed or refuted by various other investigators. 
This has been our motivation for approaching the search from a different direc- 
tion: applying the Virial Theorem to extant data. We conclude that the vertical 
density profile of the disk is not in a state of equilbrium and, therefore, that the 
Poisson-Boltzmann approach is inappropriate and it thereby leads to indefensible 
conclusions. 

Subject headings: stars: kinematics and dynamics — solar neighborhood — dark 
matter 



1. Introduction 



The mass distribution of our Galaxy can be derived from a knowledge of forces both 
parallel and perpendicular to its midplane. Our knowledge of the parallel forces comes by 
inference from measurements of the generally flat rotation profiles of other spiral galaxies 
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and by direct measurements of Milky Way stellar motions that are used to estimate galactic 
rotation. The conventional explanation for the observations is that dark matter (DM) in a 
halo somehow "conspires" with visible matter to keep the rotation profiles flat. 

Studies of forces perpendicular to the galactic midplane have also led to DM hypotheses, 
but this DM is confined to the disk where it purpor tedly has a h igher density than the halo 
DM. The search for disk DM, started 95 years ago by lOpikl ( 119151 ). h as a controversial history 
It has lasted even longer than the 75 year futile search for Vulcan ( iBaum fc Sheehanlll997l ). 



Kapteynl (119221 ) was the first to apply the kinetic theory of gases to i nterpr e t star 



counts and estimate the gravitational field of the disk in the vicinity of the Sun. lOortl (119321 ) 
refined the theory and wrote the seminal paper on the topic. The gravitational potential 
can be estimated from observations of the midplane or column densities and the velocity 
distributions of all observable stars and gases, and it can independently be estimated from 
observations of limited sets of tracer stars. A discrepancy between the two estimates may 
be indicative of DM. Numerous papers wer e published following Oort's lead, culminating 
with the disk model of iBahcall et al.l (119921 ) (BFG). The data they used to determine the 
gravitational potential comprised 13 stellar components (excluding a very small spheroid 
component) plus four gas components. They found that their data supported the hypothesis 
that disk DM exists. Yet the significance of this finding is debata ble because it has not 
been confirmed on using K giant data f rom the Hipparcos c a talogu e (jKorchagin et al 2003 ; 
Holmberg fc Flvnnl [20o3 ) . Specifically, Irlolmberg fc Flvnnl tooi ) (HF) updated the BFG 



model and proffered the most comprehensive attempt yet at modeling the galactic disk mass 
density in the solar vicinity. The HF model has 11 stellar components and the same four 
gas components as the BFG model, plus a thick disk constituent, but no disk DM. 

Another a pproach is the s tudy of the thicknes s of th e neutral hydrogen component of 
the Milky Way. Kalberlal ( 120031 ) and lKalberla et al.l ( 120071 ) analyzed flaring data in a Galaxy 
model with 6 gaseous and 4 stellar components (including a central bulge and bar). To ex- 
plain the observations, they required DM in a thick disk and i n a ring - p l us the usual DM 
halo. Using a computer simulation of the evolution of galaxies, iRead et al.l ( 120081 ) concluded 
that our Galaxy includes disk DM . Conversely, analyse s of gravitational microl ens obser- 
vations toward the galactic center ( IPopowski et a/.ll2005l ; lHamadache et all 120061 ) found no 
evidence for DM in the disk; and, alternatiy ely, the semblance of a similar, yet fictitious, DM 
disk is predicted by MOND rtMikroml 120011 ) . Disk DM, like halo DM, is a broadly accepted 
working hypothesis, but it remains just a hypothesis that demands further investigation. 
Independently being able to confirm or deny the existence of a DM disk would be important 
to fundamental physics and cosmology, so we have devised a technique to analyze the z and 
v z distributions of Milky Way stars and gases by applying the Virial Theorem (Section 12. 3p 
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to the HF disk model (Section 13.11) . 



2. Theory 

2.1. The Virial Theorem 

We adopt the BFG toy model for which the mass density of a uniformly flat Milky Way 
disk is p(z) = p(—z). In the vicinity of the Sun, at distance r = r from the Galaxy rotation 
axis, the disk mass surface density (also called the column density) is 



p(z) dz; 

J — oo 

its kinetic energy surface density is 

T = S(T 2 /2 (1) 

(which, in effect, defines the variance a 2 of the vertical velocity v z ); and its gravitational 
energy surface density is 

/oo /*oo 
pty' zdz = / pgzdz, (2) 
-oo J — oo 

where \l/ is the gravitational potential and g is the gravitational specific force (acceleration). 
For {z/r Q ) 2 <C 1, the gravitational energy between horizontal layers parallel to the midplane 
is effectively independent of the gravitational energy of the orbit about the Galaxy rotation 



axis. Then if the disk does not secularly disperse, the Virial Theorem (IGoldsteinl 11959 



Collins! 119781 ) applies: 



2T + V = 0, (3) 

where the overbar is the time averaging operator. If the vertical mass column is in a state 
of equilbrium [ p(z) = ], then the instantaneous value of IT + V is 0. If instead the mass 
column is not in a state of equilibrium, then there is an interchange of energy density back 
and forth between T and V, and consequently 2T + V oscillates about 0. 



2.2. The Poisson-Boltzmann Equation 

For a disk vertical mass column that is in an equilbrium state, the Poisson equation 
applies: 

V"(z) =4ttG P (z), (4) 
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where G is the gravitational constant. The boundary conditions are ^(0) = ^'(0) = 0. If an 
assembly of disk stars is isothermal, then a 2 is independent of z, and if it is not isothermal, 
then it is composed of multiple sub- assemblies that are isothermal. The density distribution 
of all disk stars and gas clouds, ultimately partitioned among N isothermal components, is 
the sum of iV Boltzmann distributions, 

N N 

P {Z) = Y,Pn(z) = J>n(0) exp [-*(*)/<7*]. (5) 
n=l n=l 

Combining the Poisson (T4J) and Boltzmann (jSJ) equations gives the Poisson-Boltzmann equa- 
tion, 

N 

y" = 47rGY,Pnexpi-y/a 2 n ]. (6) 

ra=l 

Spitzer derived an analytic solution for the Poisson-Boltzmann equation with a single 
(n — 1) isothermal component, 

tf' 7 = 47rGp(0) exp [-^/a 2 ] = 2k 2 a 2 exp [-^/a 2 ], (7) 

where 

k 2 o 2 = 2nGp(0). (8) 
Multiply both sides of Eq. [7] by ty', and then integrate, 

($') 2 = 4jfeV[l -exp(-^/a 2 )]. (9) 



Define 
then 

so the LHS of Eq. M is 
and the RHS of Eq.[9]is 



u- 2 = p(z)/p(0) = exp(-^/a 2 ); 
* = 2a 2 \nu, 

($') 2 = Aa\u'/u) 2 , (10) 



4JfeV 4 [l - exp (-^/a 2 )] = 4fcV[l - u~ 2 ]. (11) 
Equating Eqs. [10] and ITT] 

, = kdz. 
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Integrating this equation from u(0) = 1 to u{z) gives u(z) = cosh kz, and then the Spitzer 
solution, 

p(z) = p(0)u- 2 = p(0) sech 2 kz , (12) 

and 

V(z) = 2a 2 lnu = 2a 2 In cosh kz . (13) 

This solution conforms with the steady state Virial Theorem (2T + V = 0) because, by 
Eqs. [12 and OQ 

/oo 
sech 2 kzdz = 2p(0)/fc, (14) 
-oo 

2T = 2p(0)a 2 /k; (15) 

and, by Eqs. M, EES El and UM, 



, s n f°° ,o, rffln cosh /cz) 
1/ = -2p OW 2 / sech 2 ^^ ^cfe 

i-oo ^ 



— oo 

/■OO 

-2p(0)a 2 / sech 2 ^ tanhkzkzdz = -2p(0)a 2 /k = -2T. (16) 



We are aware of no analytic solution to Eq. |6] for N > 1, but it can be integrated 
numerically to find a unique, steady state solution. [This is not a trivial point because 
no such solutions can exist for 3-D Poisson-Boltzmann in tegrations such as modeling the 



structure of an isothermal star ( iPestana &: Eckhardtl 120071 ).] We are accordingly poised to 



examine the problem by using the Virial Theorem. The Virial Theorem approach has an 
advantage over Poisson-Boltzmann integration because it merely requires that 2T + V = 0, 
so an oscillatory solution is feasible, whereas the Poisson-Boltzmann integration requires 
that 2T + V = and so only a time invariant solution is allowed. 



2.3. The Vertical Structure of the Disk 

An infinitesimal layer of matter, p(z) dz, is the source of a specific force 

g = 2nGp(z) dz 

in the direction of the layer. Thus the gravitational energy surface density between infinites- 
imal layers of matter at z and z + a is 

— ga p(z + a) = —2nGa p(z) p(z + a) dz, 
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and the gravitational energy surface density due to all layers that are separated by distance 
a is 

/oo poo 
p(z)p(z + a) dz = —2irGY, 2 a / p(z)p(z + a) dz = —2nGT l 2 q(a)a, 
-oo J —oo 

where p(z) is the frequency function 

p(z) = p(z)/E, (17) 
and q(a) is the frequency function given by the convolution 

/oo 
p(z) p(z + a) dz. (18) 
-oo 

The total gravitational energy surface density is then 

POO 

V = -2vrGS 2 / q(a)ada. (19) 
Jo 

The lower limit of this integral is because a > to avoid double counting. 
Combining Eqs. [Tj [3] (not time averaged), [T7] and dH] gives 

p(0)o* 



POO 

2p(0) / q(a)ada. (20) 
Jo 



Since both sides of Eq. |20] are dimensionless, they need not have the same units of length. 
Then for computational convenience we implicitly define the unit of length on the RHS (and 
equations leading to it) by setting 

p(0) = 1/2, (21) 
so Eq. [20] becomes the Virial theorem criterion for a state of equilibrium, 

R = Q, (22) 

where 



p(oy 2 

and 

/■oo 

Q= / q(a)ada. (24) 

In Eq. 123], ^/p(0) and cr 2 /GE retain the dimension pc. Using Eqs. [81 and [HI this ratio 
for the Spitzer solution is 

R s = 1/2. (25) 
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We have examined various trial frequency functions, p(z), using them in Eqs. [TBI and 
[Ml to evaluate Q. These functions are not arbitrary. Each must be bell-shaped: an even 
function of z (that is, p(z) = p(—z)), with a single peak at p(0) = 1/2 (that is, zp'(z) < 0) 
and, of course, p(z)dz = 1. Moreover, we assume that p(z) and all its derivatives are 
continuous, and that p(z) is not a pathological frequency function, meaning that its n-th 
moment J^° oo z n p(z) dz exists. Therefore, for large z, p(z) ~ exp(— az 2m ), where a is a 
positive constant and m > 1 is an integer. With these provisos, it turns out that for all 
plausible frequency functions, 

Q « 1/2. (26) 
We offer two simple examples (for which m — 1): 

• For the Spitzer function (see Eqs. [121 ITTlandl2~l| with k = 1), 

Ps{z) = - sech 2 z, (27) 

which applies exactly for a single component model, but still might be a suitable 
approximation for a multiple component model. For the Spitzer function, 

Qs = 1/2, (28) 
so, as expected (see Eq. [25]) , Rs = Qs- 

• For the Gaussian frequency function, 

(29) 
for which 

Q G = y/2/n = 0.45. (30) 

These are merely examples. The exact value for Q is not important for this study. Eq. [221 is 
entirely adequate to restate the condition for equilibrium, Eq. [22J as 

R = Q& 1/2. (31) 



Pg{z) 



2 6XP 



■ z \ 2 

- Tl .2, 



3. Disk Mass Models 
3.1. The HF Model 



The 15 components of the HF disk mass model were itemized by iFlynn et all 12006 
From these, we calculate 



15 



E = J^^i = 49.3 pcT 2 M , (32) 
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ttGS 2 = tt x 49.3 2 pc" 4 (GM Q ) M q = 34.4 x 1(T 27 pc'V 2 M , (33) 

14 

a 2 = ^(lysy, 2 = 587 km 2 s" 2 = 614 x l(r 27 pc 2 s- 2 , (34) 
l 

15 

P(O) = = 0.0914 pc- 3 M , (35) 

i 

p(0)a 2 = 56.1 x 1(T 27 pc -1 s -2 M , (36) 

and (see Eq. |23I) . 

Rhf = 1-63, (37) 

which is too large compared with Q « 1/2 for an equilibrium solution to exist. Hence, a 
Poisson-Boltzmann integration using the HF model, or any other such model, is an inap- 
propriate and futile exercise. Nevertheless, HF performed this integration without disk DM 
and found the result to be compatible with vertical density profiles of K giant tracers. We 
look at the same data and conclude that the HF model can be compatible with a Poisson- 
Boltzmann steady state solution only if disk DM exists. For example, suppose that DM has 
the same vertical profile and v z variance as visible matter, but that visible matter accounts 
for only 30% of the total surface mass density. Then p(0) and S increase by the factor 10/3, 
and thus Rhf is modified to 1.63 x 3/10 = 0.49, which is a plausible value for Q. How- 
ever, the revised model no longer satisfies the Poisson-Boltzmann equation and the change 
increases the vertical accelerations of the tracer stars by the factor 10/3 and decreases their 
scale height by 70%, and such a model conflicts drastically with K giant vertical profiles. 
Introducing DM solves one problem but creates another. The overall problem can only be 
resolved by abandoning the requirement that the disk z profile is in a steady state or that 
gravity is Newtonian. 



3.2. The BFG Models 

BFG conceived and examined ten disk models, nine of them containing varying formu- 
lations with DM. Their Table 4 lists their observed model in Row 1 [pi(0) = 0.1026 pc~ 3 M 
and Si = 49.8pc~ 2 M ], and their "best fit" model (observed plus dark matter) in Row 3 
[p 3 (0) = 0.2596 pcr 3 M and £ 3 = 83.9 ptr 2 M ]. For the Row 1 model, the ratio of Pi(0)/£ 2 
to p HF (0)/T, 2 HF is (0.1026 x 49.3 2 )/(0.0914 x 49.8 2 ) = 1.100, and for the Row 3 model, the 
ratio of p 3 (0)/S 2 to p HF (®)/Y? HF is (0.2596 x 49.3 2 )/(0.0914 x 83.9 2 ) = 0.981. If a 2 does not 
change between the HF and BFG models, R\ = 1.79 and R3 = 1.60. Like Rhf = 1-63, R\ 
and R3 are incompatible with R = Q w 1/2. 
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4. An Oscillating Disk 

A Poisson-Boltzmann solution is feasible only if 2T + V — 0, but this requirement is 
incompatible with observations; compelling a steady state solution is tantamount to com- 
pelling the existence of DM or non- Newtonian gravity. Suppose instead that a disk column's 
vertical profile varies periodically with time. Then 2T + V oscillates about 2T + V = 0, so 
that at any epoch, 2T + V could be substantially greater or less than zero. In other words, 
R oscillates about Q. We use this broadened aspect to interpret the HF data, and uncover 
interesting possible ramifications. 

The analysis so far has been quantitative, but from here on it will be descriptive and 
speculative - and accordingly less defensible. Using dimensional analysis, we reckon that the 
z oscillation period is of the order [Gp(O)] -1 / 2 = 50 million years. The estimated period 
is roughly a quarter of a galactic year, but this estimate is inexact because it is based 
on dimensional considerations only. It would have been larger if we had used the mean 
midplane density p(0) over one cycle instead of the current midplane density (which we 
believe to be near the maximum value), and there undoubtedly are dimensionless factors 
that would arise in a precise analysis, decreasing the estimated period by as much as a factor 
of 10, but increasing it no more than fourfold (each disk star orbiting the galactic center has 
a retrograde precession caused by the other disk stars, so the z oscillation period must be 
less than the orbital period). What we find puzzling about the stellar orbits is the cause of 
the symmetric clustering of their ascending and descending nodes. 

Disk column oscillations are exchanges between gravitational, V oc — £ 3 /p(0), and ki- 
netic, T oc Ha 2 , energy surface densities. Since £ does not vary with time, oscillations in 
p(0) entail oscillations in a 2 ; and the p(0) and a 2 variations are in phase, so that the peaks 
in midplane mass and kinetic energy densities coincide. Eq. [37] tells us that at the current 
oscillation phase p(0)<r 2 is about three times higher than its mean value, but it does not tell 
us how close p(0)<r 2 is to its peak value. At the peak, there would be no net flow of stars 
toward or away from the midplane. Unaware of any evidence for such a flow, w e suppose that 



p(0)cr 2 is close to its peak, and corroborate this with data from Table A.l of iFamaey et al. 



( 120051 ). selecting single giants (* tag in Column 24) in the smooth background (g — 3 in 
Column 34) that are in the two zones defined by 100 < \z\ < 150 pc. There are 467 stars in 
the northern zone and 224 stars in the southern zone. We calculate that the two zones are 
approaching each other at 1.3 ± 1.4 km s -1 . The Virial Theorem requires an oscillation with 
substantial energy flow, but these numbers tell us that there is negligible net mass flow now 
(at least for the giants), so the present epoch appears to be one of maximum p(0)<r 2 . With 
little net mass flow toward or away from the midplane, the v z velocity distributions may 
appear to be isothermal, but at epochs of maximum flow, the net flow velocities would be of 
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the order S[G/p(0)] 1 / 2 = lOkms 1 . Considering its imprecision, this is not much different 
from guf — 24kms _1 . A net outward flow would be obvious in about 10 million years. 

We extend our z coordinate system to a standard cylindrical coordinate (r, z) system. 

1/2 

Then we expect oscillations elsewhere on the disk with periods proportional to p(r, t?, 0) 
Because these mean midplane densities vary, so do the z oscillation periods. Just small 
differences in p(r, 0) between neighboring disk columns can effect large changes in their 
relative z oscillation phases and, therefore, large differences in their instantaneous midplane 
densities p(r, 0). 

The increased midplane mass and kinetic energy densities caused by oscillations are 
conducive to the formation of stars, including the massive O and B stars whose luminosity 
lifetimes would be relatively ephemeral compared with the disk oscillation periods. We 
therefore conjecture that the disk spirals are where the midplane mass and kinetic energy 
densities are now high, so stars are being formed there at relatively high rates. Vertical 
oscillations also occur in the gaps between the spirals, but their midplane mass and kinetic 
energy densities are now low (and their column gravitational energies are high) because of 
phase differences, so brighter stars that were created when their midplane mass and kinetic 
energy densities had been high have sinced waned. A computer simulation might tell us 
the shape that p(z) takes in the gaps, for we do not even know whether the frequency 
function remains unimodal. But neither do we know whether our conjecture concerning the 
formation of spirals is correct; nor, indeed, can we be sure that the disk actually oscillates, 
for conceivably the HF data can be explained with a non-Newtonian gravitational theory. 
However, we maintain that the conventional Poisson-Boltzmann interpretation of the data is 
demonstrably wrong, and we offer the disk oscillation hypothesis as one plausible alternative. 



5. Conclusion and Implications 



We have used the Virial Theorem to study the problem of determining the density of 
matter near the Sun because it is less restrictive than the Poisson-Boltzmann approach. In 
fact, we have shown that all the comparisons of tracer star de nsities with Pois so n-Bol t zman n 
models or with equ ili brium models in general, for exampl e iKaptevnl (119221) ; lOortl (119321 ); 
Bahcall et all f)l992h : iKalberlal (ho03h : iHolmberg fc Flvnnl (l200ol . booi ). have been wrong. 
The reason is that the disk is not in a state of equilibrium (2T + V ^ 0). Furthermore, any 
model made compatible with 2T + V = by adding DM will then be incompatible with 
Poisson-Boltzmann integration models, and with observed tracer star densities. No evidence 
can be f o und f or DM by using this approach. Explaining the evidence (e.g., the BFG, 
Kalberlal (120031 ) and HF disk models, plus tracer star densities) requires disk oscillations or, 
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perhaps, non-Newtonian gravity. A cursory inquiry into the possibility of disk oscillations 
without DM leads to intriguing hints concerning the structure and dynamics of a galactic 
disk. 

We cont end that any as tr ophysical argument pre d icated on the existen c e of a DM disk 
for e x ample IStotherd (119981); ISanchez-Salcedol (19991) ; Ide Boer et al.l (120051 ) ; iKalberla et al. 



( 2007 ); Revaz et al. ( 20091 ); Gates fe Gyuk ( 2004 ) - does not warrant consideration. We reject 



the recently proposed model for the f ormation of disk galaxies within the A CDM framework 
because it would produce a DM disk ( Read et al. 2008 ; Purcell et al. 20091) . Finally, there is 



no fic titious semblance of a DM disk as predicted by MOND (lMilgromll200ll ; iBienayme et al. 
20091 ). so our analysis argues against MOND. 



We thank the referee for positive comments and constructive suggestions, and especially 
for drawing out attention to a key numerical error and encouraging us to resolve its effect. 
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